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ABSTRACT 

We discuss experimental data on particle yields and particle spectra obtained in heavy ion col- 
lisions in a very broad energy range from SIS/GSl {^/s ~ 2 GeV) through AGS/BNL {^/s ~ 5 
GeV) up to SPS/CERN (v^ ~ 20 GeV) and RHIC/BNL (^/i ~ 130) GeV. We argue that in 
this broad energy range hadronic yields and their ratios resemble a thermal equilibrium pop- 
ulation along a unified freeze-out curve determined by the condition of fixed energy/particle 
~ 1 GeV. At RHIC and top SPS, thermal parameters are consistent within error with the 
critical conditions required for deconfinement. This, together with the particular distribution 
of strangeness within a collision fireball, could indicate that chemical equilibrium is a direct 
consequence of parton to hadron transition, which populates a state of maximum entropy. At 
lower energies equilibration in A-A collisions should appear through hadronic interactions and 
rescatterings. 
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1 Introduction 



The ultimate goal of ultrarelativistic nucleus-nucleus collisions is to study the properties of 
strongly interacting matter under extreme conditions of high energy density Quantum 
Chromodynamics (QCD) predicts that strongly interacting matter undergoes a phase transition 
from a state of hadronic constituents to a plasma of unbounded quarks and gluons (QGP) p|, |[. 
By colliding heavy ions at ultrarelativistic energies, one expects to create hadronic matter under 
that conditions are sufficient for deconfinement 0, ^, ^ |, Thus, of particular relevance was 
finding experimental probes to check if the produced medium in its early stage was indeed in 
the QGP phase. Different probes have been theoretically proposed and studied in terms of 
SPS/CERN and most recently RHIC/BNL experiments. The most promising signals of QGP 
were related with particular properties of photons ^, dileptons [|], |TD[ and hadron spectra 



0, ^, 0. The photon rate was expected to be enhanced if the QGP were formed in the initial 



state. The invariant mass distribution of dileptons should be modified by in-medium effects 
related with chiral symmetry restoration ||10| , |12[. The suppression of charmonium production 
was argued to be a consequence of the collective effects in a deconfined medium [||, [T^ . 

Hadron multiplicities and their correlations are observables which can provide information 
on the nature, composition, and size of the medium from which they are originating. Of 
particular interest is the extent to which the measured particle yields are showing equilibration. 
The appearance of the QGP, that is a partonic medium being at (or close to) local thermal 
equilibrium, and its subsequent hadronization during phase transition should in general drive 
hadronic constituents towards chemical equihbrium [H, ^ H, [1^. Consequently, a high level of 



chemical saturation, particularly for strange particles could be related with the deconfined 
phase created at the early stage of heavy ion collisions. 

The level of equilibrium of secondaries in heavy ion collisions was tested by analysing the 
particle abundances p|, 10, |T6|, [13, or their momentum spectra |Tl|, EH, ETl. In the first 



case one establishes the chemical composition of the system, while in the second case additional 
information on dynamical evolution and collective fiow can be extracted. 

In this article we will discuss and analyse the experimental data on hadronic abundances 
obtained in ultrarelativistic heavy ion collisions, in a very broad energy range starting from 
RHIC/BNL (v^ = 130 A GeV), SPS/CERN (^i ^ 20 A GeV) up to AGS/BNL (v^ ~ 5 
A GeV) and SIS/GSI (a/s ~ 2 A GeV) to test equilibration. We argue that the statistical 
approach provides a very satisfactory description of experimental data covering a wide energy 
range from SIS up to RHIC. We discuss the unified description of particle chemical freeze-out 
and the excitation function of different particle species. Introducing, in addition to thermal, 
also the transverse collective motion, the systematics of thermal freeze-out is also presented. 



2 Initial conditions in A— A collisions and deconfinement 

In ultrarelativistic heavy ion collisions, the knowledge of the critical energy density ec required 
for deconfinement is of particular importance as well as the equation of state (EoS) of strongly 
interacting matter. The value of ec is needed to establish the necessary initial conditions in 
heavy ion collisions to possibly create the QGP, whereas EoS is required as an input to describe 
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the space-time evolution of the colhsion firebalL 

Both of these these pieces of information can be obtained today from first principal calcu- 
lations by formulating QCD on the lattice and performing numerical Monte-Carlo simulations. 
In Fig. 1 we show the most recent results of lattice gauge theory (LGT) for energy density 



and pressure p3|- These results have been obtained in LGT for different numbers of dynam- 
ical fermions. The energy density is seen in Fig. 1 to exhibit a typical behaviour in a system 
with a phase transition: an abrupt change in the very narrow temperature range.Q The corre- 
sponding pressure shows a smooth change with temperature. In the region below the basic 
constituents of QCD, quarks and gluons, are confined within their hadrons and here the EoS 
is well parametrized by the hadron resonance gas. Above Tc the system appears in the QGP 
phase where quarks and gluons can penetrate distances that substantially exceed a typical size 
of hadrons. The most recent results of improved perturbative expansion of thermodynamical 
potential in the continuum QCD are showing that at some distance above Tc the EoS of 
QGP can be well described by a gas of massive quasi-particles whose mass is temperature de- 
pendent. In the near vicinity of Tc the relevant degrees of freedom were argued to be described 
by the Polyakov loops ||2^ . 

Lattice Gauge Theory predicts that in two-flavour QCD the critical temperature Tc = 
173 ± 8MeV and the corresponding critical energy density ec = 0.6 ±0.3 GeV/fm^ are required 
for chiral phase transition . The value of ec is relatively low and quantitatively corresponds 
to the energy density inside the nucleon. 




Figure 1: Pressure P and energy density normalized to temperature in fourth power, versus 
temperature normalized to its critical value. The calculations were done within LGT for 
different numbers of flavours [^]. The values of the corresponding ideal gas results are indicated 
by arrows. 

The initial energy density reached in heavy ion collisions can be estimated within the Bjorken 



model 26|. From the rapidity distribution of protons and their transverse energy Et measured 



in nucleus-nucleus collisions the initial energy density eo is determined from 



^We have to point out, however, that in the strictly statistical physics sense, a phase transition can only 
appear in the limit of massless quarks. 
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TxK'^ tq ay 

where the initially produced collision fireball is considered as a cylinder of length Tody and 
transverse size R ~ A^/^ . Inserting for txR^ the overlap area of colliding Pb nuclei together 
with initial time an tq ~ 1 fm, and using an average transverse energy at midrapidity measured 
at the SPS (v^ = 17.3 GeV) to be 400 GeV |3, one obtains 

el^^{TQ ~ 1 fm) ~ 3.5 ± 0.5 GeV/fml (2) 

Increasing the collision energy to = 130 A GeV for Au-Au at RHIC and keeping 
the same initial thermalization time as at the SPS, one would expect an increase of eo by 
only 50-60 %. However, at RHIC the thermalization is argued, within saturation models, to 
appear at much shorter time. The basic concept of saturation models is a conjecture that 
there is some transverse momentum scale p^at where the gluon and quark phase space density 
saturates |2^. For isentropic expansion of the collision fireball, one can relate the transverse 
energy at psat with the one measured in nucleus-nucleus collisions in the final state. The 
saturation scale also fixes the time scale tq = 1/psat- Taking the value of psat predicted in 
g|] for RHIC Psat — 1-13 GeV, one gets tq — 0.2 fm and the corresponding initial energy 
density eo — 98 GeV/fm^. The estimate of eo, however, strongly depends on the value of psat, 
which is model-dependent. The McLerran-Vanugopalan model ||3^, for instance, predicts the 



value of e^^^^ ~ 20 GeV/fm^, which agrees with the predictions of |31|. At the SPS the 
saturation model described in ^9] leads to e^^^ ~ 16 GeV/fm^, a much higher value than 



given in Eq. (2). It is thus clear that there are large uncertainties on the value of the initial 
energy density reached in ultrarelativistic heavy ion collisions. In Pb-Pb collisions at the SPS 
according to the models one gets 2.5 GeV/fm^ < e^^^ < 16 GeV/fm^ whereas in Au-Au 
colhsions at RHIC, 20 GeV/fm^ < e^^^ < 100 GeV/fm^. 

The dominant component of the partonic medium produced in ultrarelativistic heavy ion 
collisions at RHIC and even at the SPS is gluons. The energy density of gluons in thermal 
equilibrium scales with the fourth power of the temperature e = gT^, where g denotes the 
number of degrees of freedom. For an ideal gas, g = 167r^/30; in an interacting system, the 
effective number of degrees of freedom g is smaller. The results of LGT shown in Fig. 1 indicate 
deviations from the Boltzmann limit by 20-25 %. Relating the thermal energy density with 
the initial energy density discussed above, one can make an estimate of the initial temperature 
reached in heavy ion collisions. For the SPS this gives a temperature in the range 200 MeV < 
T < 330 MeV and at RHIC 400 MeV < T < 600 MeV. 

Comparing the initial energy density expected in heavy ion collisions with LGT results, it 
is clear that the initial energy density, at RHIC, exceeds by far the critical value. Thus, the 
necessary conditions to create the partonic medium in a deconfined phase are reached at RHIC 
as well as at the top SPS. Large energy density is, however, still not sufficient to create a QGP. 
The distribution of initially produced gluons is very far from being thermal, thus the system 
needs enough time to equilibrate. Recently, it was rigorously shown in the framework of 
perturbative QCD and kinetic theory that the equilibration of partons happens indeed at the 
LHC and most likely at RHIC. Previous microscopic study within the Parton Cascade Model 
has also suggested that thermalization can be reached at lower SPS energy |3^. Here, however. 
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Figure 2: Comparison of the experimental data on different particle multiplicity ratios obtained 
at RHIC at {^/s — 130) GeV with thermal model calculations for T — 175 MeV and /jLb — 51 
MeV. 

it is not clear if models inspired by the perturbative QCD are indeed applicable at this relatively 
low collision energy. 

Admitting QGP formation in the initial state in heavy ion collisions one could expect that 
the thermal nature of the partonic medium could be preserved during hadronization. Conse- 
quently, the particle yields measured in the final state should resemble the thermal equilibrium 
population. In the following, we present the most recent results related with the question of 
equilibration of secondaries in heavy ion collisions and discuss its possible relation with decon- 
finement. 



3 Statistical model and particle multiplicity 

The basic quantity in the statistical model description of thermal properties of hadronic matter 
is the partition function Z{T, V). In the Grand Canonical (GC) ensemble, 

Z^^{T, V, iiq) = Tr[e-^(^-Si'^«i'2*)], (3) 

where H is the hamiltonian of the system, Qi are the conserved charges and /xg. is the chemical 
potentials that guarantees that the charge Qi is conserved on the average in the whole system. 
Finally /? = 1/T is the inverse temperature. 

In the strongly interacting medium, one includes the conservation of electric charge, baryon 
number and strangeness. Thus, the partition function depends in general on five parameters. 
However, only three are independent, since the isospin asymmetry in the initial state fixes the 
charge chemical potential and strangeness neutrality conditions eliminate the strange chemical 
potential. On the level of particle multiplicity ratios derived from the partition function, we 
are left with only temperature T and baryon chemical potential /x^ as independent parameters. 
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Figure 3: On the left, the ratio of kaon to pion measured in Au-Au colhsions at 1 A GeV 



IJ]; the broken hne represents the statistical model results The rights-hand figure shows 
statistical model predictions for yield/participants in A-A collisions at 40 A GeV normalized 
to the corresponding value in pp collisions. 



The Hamiltonian is usually described by the hadron resonance gas, which contains the 
contributions from all mesons with masses below 1.6 GeV and baryons with masses below 2 
GeV. In this mass range the hadronic spectrum is well established and the decay properties of 
resonances are known. This mass cut in the contribution to partition function limits, however, 
the maximal temperature to Tmax < 190 MeV, up to which the model predictions could be 
trustworthy. For higher temperatures the contributions of heavier resonances are not negligible. 

In the high density regime the repulsive interactions of hadrons should also be included in 
the partition function. To incorporate the repulsion on short distances between hadrons one 
usually uses a hard core description by implementing the excluded volume corrections. In the 
thermodynamically consistent approach [34| these corrections lead to a shift of baryon chemical 
potential. Finally, the widths of resonances and their decay into lighter particles have to be 
included in the statistical model when calculating particle multiplicities [36[] . 

The statistical model, described above, was applied to Pb-Pb collisions at top SPS energy 
1^, [18|, 0. The model was compared with almost all experimental data obtained by NA44, 



NA49 and WA97 Collaboration. Hadron multiplicities ranging from pion to omega and their 
ratios were used to verify if there is a set of thermal parameters (T, /x^) that simultaneously 
reproduces all measured yields. A detailed analysis has shown that choosing a temperature 
T = 168 ± 4 MeV and a baryon chemical potential hb = 266 ± 8 MeV, the statistical model 
with only two parameters can indeed describe seventeen different particle multiplicity ratios 
within an accuracy of one to two standard deviations. One could thus conclude that, with 
respect to the statistical operator formulated for equilibrium hadron resonance gas, the experi- 
mental data at the SPS are showing a high level of chemical equilibration. The natural question 
arising here would be to what extent this statistical operator provides a unique description of 
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the data. This question was addressed in the hterature and two distinct models have been 
examined ^ |3^. In the authors analysed the possible influence of in-medium effects on 
the chemical equilibrium description of particle yields at the SPS. In [0 the non-equilibrium 
scenario of explosive hadronization of a QGP fireball was proposed. Both these models and 
particularly are showing satisfactory agreement with SPS data, however with larger devia- 
tions for multistrange particles. In our discussion we concentrate on equilibrium description of 
particle production since only this approach, as will be demonstrated, provides the systematic 
agreement with almost all heavy ion data from SIS to RHIC. 

The chemical freeze-out temperature, found from a thermal analysis ||18|, |19[ of experimental 
data in Pb-Pb collisions at the SPS is remarkably consistent, within errors, with the critical 
temperature Tc ^ 173 ± 8 MeV obtained from lattice Monte-Carlo simulations of QCD at 
vanishing baryon density p3|. Thus, the observed hadrons seem to be originating from a 
deconfined medium and the chemical composition of the system is most likely to be established 
during hadronization ^, |^ . The observed coincidence of chemical and critical conditions in 
the QCD medium, if indeed valid, should be seen also in heavy ion collisions at higher collision 
energies, in particular at RHIC. 
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Figure 4: Particle yields per participant in Pb+Pb relative to p+Be and p+Pb collisions 
centrality dependence. The data are from WA97 HO] and NA57 H^] Collaborations. 



The equilibrium statistical model was recently applied to Au-Au collisions at a/s = 130 GeV 
at RHIC The results of the STAR, PHENIX, PHOBOS and BRAHMS Collaboration for 



different particle multiplicity ratios have been used to test chemical equilibration at RHIC. In 
Fig. 2 we show the comparison of the thermal model with experimental data. One sees that 
the overall agreement is very good. Most of the data are reproduced by the model within the 
experimental errors. The largest deviations are seen in the ratios of K*^ jhr and K^^/h' but 
they are still on the level of one standard deviation.^ In Au-Au collisions at \/s= 130 GeV 
the chemical freeze-out appears at T = 175 ± 7 MeV and /^s = 51 ± 6 MeV. The resulting 
temperature is only slightly higher than that previously found at the SPS where for Pb-Pb 

^ In the model calculations, the h~ was considered as the total number of negatively charged particles in the 
thermal fireball 
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collisions T = 168 ± 5MeV. This relatively moderate increase of temperature could be expected 
since, in the limit of vanishing baryon density, the temperature should not exceed the critical 
value required for deconfinement. The substantial decrease of the baryon chemical potential 
from fiB — 270 MeV at the SPS to fiB — 50 MeV at RHIC shows that, at midrapidity, we are 
dealing with a low net baryon density medium. 

The results for particle yields and their ratios at the SPS and RHIC shows the statistical 
order. Chemical equilibration of secondaries after hadronization is rather excluded by kinetics 



1^, ^]. Thus, the equilibrium population of hadrons would be most likely to appear since 



it was pre-established in the QGP phase. In the following we argue, however, that equili- 
bration of secondaries is not a unique signal for deconfinement, as it is also there at lower 
energies, where the initial conditions exclude QGP formation. To test equilibration in low 
energy nucleus-nucleus collisions, one needs, however, to change the statistical operator from a 
GC to a canonical C ensemble with respect to strangeness conservation. 



4 Equilibrium limit of rarely produced particles 

The conservation of quantum numbers related with U(l) internal symmetry in statistical models 
can be described in the GC ensemble only if the number of produced particles per event carrying 
corresponding quantum number is much larger than 1. In the opposite limit of rare particle 
production [^, U(l) charge conservation must be implemented locally on an event-by- 
event basis, i.e., a canonical C ensemble of conservation laws must be used. The C ensemble 
is relevant in the statistical description of particle production in low energy heavy ion 



or 



high energy hadron-hadron or e^e reactions PS| as well as in peripheral heavy ion collisions 
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The exact conservation of quantum numbers, that is the canonical approach, is known to 
severely reduce the thermal phase-space available for particle production [|^. Consequently, 
the chemical equilibrium limit of rarely produced particles is changed and it is different from 
the one obtained in the asymptotic GC limit. In order to illustrate the above change, let us 
consider the kinetics for the time evolution and equilibration of rarely produced particles by 
considering a simple example of K^K~ pair production and equilibration in the environment 
of thermal pions contained in volume V at temperature T. The production of kaons is due to 
the binary process tt+tt^ K^K^ . 



In the standard formulation [Q, the rate equation for this binary process is described by 
the following population equation: 

^ = ^(iV.^)(iV.-) - ^{Nk^){Nk-), (4) 

where G = {(Jgv) and L = {(Tlv) give the momentum-averaged cross sections for the gain 
TT+TT^ K^K~ and the loss K^K^ ti^ti^ process, respectively, and (Nk) represents the 
total number of produced kaons. 

In the above equation it is explicitly assumed that and K" are uncorrelated. To 
include possible correlations between the production of and K~ ^T|], let us define Ptj as 



the probability to find i number of and j number of K in an event. We also denote by 
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Thermal Model 



Pi the probability to find i number of K in an event. The average number of K per event is 
defined as: (A^^) = EZo^Pi- 

We can now write the following general rate equation for the average kaon multiplicities: 

'-^ = ^{N..){K-)-^j:^jP.,r (5) 

Owing to the local conservation of quantum numbers, we have: 

P = P 5- 

E^JA. = Y.^'P^^{N') = {Nf + {m'), (6) 

where {6N'^) represents the event-by-event fluctuation of the number of K^K~ pairs. Note 
that we always consider abundant vr"*" and vr" so that we can neglect the number fluctuation of 
these particles and the change of their multiplicities due to the considered processes. 

Following Eqs. (^) and the general rate equation for the average number of K^K^ pairs 
can be written as: 

^^ = f{iV..){iV.-)-i(iVi). (7) 

For abundant production of K^K~ pairs where {Nk) ^ 1, [N]^) ~ (Nk)"^, and Eq. (|^ 
obviously reduces to the standard form: 

^ ^{N..){N^-) - ^{Nk^){Nk-). (8) 
8 
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Figure 5: Compilation of chemical freeze-out parameters from SIS to RHIC. The broken line 
represents a phenomenological condition of chemical freeze-out of fixed energy /particle ~ 1 
GeV[0. 

However, for rare production of K^K~ pairs where {Nk) ^1, the rate equations and 
are no longer valid. We have instead {Nj^) {N^), which reduces Eq. (0) to the form 

^«f{iV..>(iV.-)-^{iV,). (9) 

Thus, the limit where (Nk) -C 1, the absorption term depends on the pair number only linearly, 
instead of quadratically for the limit of (Nk) ^1. It is thus clear, that the time evolutions and 
equilibrium values for kaon multiplicities are obviously different in the above limiting situations. 

In the limit of large (Nk) the equilibrium value for the number of K~^K~ pairs, which 
coincides with the multiplicity of and K~, is obtained from Eq. (H) as. 



(^i^)e'f = ^mlTK,iMK/T), (10) 



it is thus described by the GC result with vanishing chemical potential, in our example, because 
of the strangeness neutrality condition. 

In the opposite limit, where (Nk) -C 1, the time evolution of pion multiplicity is described 
by Eq. (H), which has the following equilibrium solution: 



(Nk) 



c 

eq 



^ MI+TK2{Mk+/T) 



2tt 



^ MI.TK2{Mk-/T) 



2n 



The above equation demonstrates the locality of strangeness conservation. With each 
the K~ is produced in the same event in order to conserve strangeness exactly and locally. 
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Figure 6: The left-hand figure shows the total number of pions per participant ((vr)/A^part)5 47r 
data in A-A and p-p collisions versus energy. Data at lower energies in A-A as well as in p-p 
collisions are from The RHIC results are from The short-dashed and dashed lines 
represent the fit to the data. The right-hand figure shows the centrality dependence of the 
K~^/K~ ratio. Data are from the STAR, NA49, E866, and KaoS Collaborations. The broken 
lines are statistical model results. 



Comparing Eqs. (|ToD and (|TT]), we first find that, for (Nk) ~ 1, the equilibrium value is far 
smaller than what is expected in the opposite limit. We also note that the volume dependence 
differs in the two cases . The particle density in the (Nk) ^ 1 limit is independent from V 
whereas in the opposite limit the density scales linearly with V^.^ 

The results of Eqs. (|T0|) and ( ]TT| ) correspond to two limiting cases of asymptotically large 
and small kaon multiplicity in heavy ion events. In order to find the equilibrium solution valid 
for any arbitrary number of kaons per events one needs to formulate a kinetic equation for the 
probability instead of the multiplicity of produced particles [ p]] . 

Let Pn{t) (0 < n < 00) be the probability function for the production of n, K^K~ pairs. 
The probability P„ tends to increase in time owing to the transition from n — 1 and n + 1 states 
to n. It also tends to decrease since the state n makes a transition to + 1 and n — 1. With the 
transition probability per unit time from n ^ n + 1 given by GNt^+Nt^-V~^ and from n — l^n 
described as LV~^, one can formulate the general iterative master equation for the probability 



function as 41 



dPn 

dr 



V ' 
L 



-N^- (-Pn-l ~ Pn 



V 



n'Pn -{n + lyPn 



+1 



(12) 



The above equation is equivalent to the general rate equation described by Eq. (J^). However, 

^ In the application of the statistical model to particle production in heavy ion collisions, the volume of 
the fireball scales with the number of participating nucleons iVpart- Thus, in terms of the canonical model, 
depending on the abundances of produced U(l) charged particles, one obviously expects different centrality 
dependences. 
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contrary to Eq. it can be solved for the equilibrium limit corresponding to an arbitrary 
number of produced kaons. Indeed, converting first the above equations for P^s into a partial 
differential equation for the generating function g{x, r) = J2'^=o x"' Pn{T) one obtains 

^^^ = ^il-x){xg" + g'-eg), (13) 

where g' = dg/dx and ^/e = (Nk)^^ given by Eq. (|TUp. The equilibrium solution for gcqlx) is 
obtained as 



9eq{x) = JJ^-J^Iq{2^), (14) 

where the normalization is fixed by g{l) = I]-Pn = 1- 

The equilibrium value for the probability function is expressed by 

^"'^^ = 7o(V^)W ^^^^ 

which converts to a Poisson distribution only in the limit of the large argument of the Bessel 
function. 

The result for the average number of kaon pairs K^K~ in equilibrium is obtained from g'{l) 
and reads: 

27r^ Io[2^mj^TK2{MK /T}\ 

The above equation is a general equilibrium solution, which is valid for an arbitrary value 
of (Nk) and obviously reproduces the asymptotic results described by Eqs. (p!OD and (pT]). 
This can be seen in the most transparent way when comparing two limiting situations: the 
large-and-small x (where x is the argument of the Bessel function) limit of the above equation. 

The equilibrium density corresponding to a large limit and described by Eq. ( p!0|) is a 
standard result for the particle density that can be obtained from the GC partition function 
introduced in Eq. (3). The general results described by Eq. (16) obviously require a different 
definition of the partition function, which takes into account the exact conservation of quantum 
numbers. This is the canonical partition function with respect to the charge conservation. In 
the C approach there is no more chemical potential under the trace as in Eq. (3) but, instead, 
the partition function is calculated by summing only those states that are carrying exactly the 
quantum number Q, that is 

Z§iT,V)^TrQ[e-^^]. (17) 

Following Eq. (|^) it is clear that C and GC limits are essentially determined by the size of 
(SN"^), the event-by-event fluctuations of the number of particles carrying U(l) charge. The 
grand canonical results correspond to small fluctuations, i.e. {6N'^)/{N)'^ -C 1, while the 
canonical description is necessary in the opposite limit. 
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Figure 7: Ratio of kaon to pion measured in heavy ion collision at different collisions energies. 
The left-hand figure corresponds to midrapidity data, whereas the right one represents the 
ratios, of fully integrated results. Data at SIS, AGS, SPS and RHIC are taken from |]SD|, [^ . 
The short-dashed line describes the results calculated following the freeze-out curve shown in 
Fig. 5. The parameterization of the p-p data from [E31 is indicated by the fuUs-line. 



The major difference between the canonical and the grand canonical treatment of the con- 
servation laws appears through a different volume dependence of particle densities as well as a 
strong suppression of the thermal particle phase space in the former. The relevant parameter 
to measure the suppression of particle multiplicities from their grand canonical result is seen, 
in Eq. (16), to be determined by the ratio of the Bessel functions /i(x)//o(x). For multistrange 
particles this suppression factor has a more complicated structure [^, |45[1 . 

In nucleus-nucleus collisions the absolute values of baryon number, electric charge and 
strangeness are fixed by the initial conditions. Modelling particle production in statistical ther- 
modynamics would in general require the canonical formulation of all these quantum numbers. 
A detailed analysis |4^, however, has shown that in heavy ion collisions only strangeness 
should be treated exactly, whereas the conservation of baryon and electric charges can be de- 
scribed by the appropriate chemical potentials in the grand canonical ensemble. 

In heavy ion collisions the number of produced strange particles depends on the collision 
energy and centrality of these collisions. At low collision energies, at SIS/GSI for example, 
the average number of strange particles produced in an event is much smaller than 1. Thus, 
here we are in the asymptotic regime of canonical ensemble. Figure 3a shows the experimental 
data on yield per participant Apart as a function of Apart measured in Au-Au collisions at 
Eiah ~ 1 A/GeV 1^^. The data are compared with the results of the canonical statistical model 
shown by the dashed-line. The thermal parameters, the temperature and the baryon chemical 
potential were chosen in such a way as to reproduce measured particle multiplicity ratios of 



strangeness neutral particles The volume parameter in the statistical operator is assumed 



to scale with the number of participants. The results in Fig. 3a clearly indicate that both the 
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Figure 8: Particle ratios in A-A versus energy. Data at the SPS are fully integrated NA49 
results (left-hand figure). The right-hand figure shows WA97 data at midrapidity. The cor- 
responding ratio at the top AGS was obtained from E810 results on measured in Si-Pb 
collisions in the rapidity interval 1.4 < y < 2.9 normalized to the full phase space values 



of 7r+ and K yield obtained in Si-Au collisions by E802 |5^. The lines represent statistical 
model results along the unified freeze-out curve from Fig. 5. 



magnitude of the yield and the strong, almost quadratic, dependence of the kaon yield on the 
number of participants is well reproduced by the canonical model. 

The importance of the canonical treatment of strangeness conservation has been shown also 
at higher collision energies, e.g., at the SPS or even RHIC, when considering the centrality 
dependence of multistrange baryons In very peripheral collisions the yield of strange 

particles is so small that the canonical description should be applied there as well. The canonical 
suppression of the thermal particle phase-space was found to increase with the strangeness 
content of the particle. 

Figure 3b shows the multiplicity /participant of Q, S, and A relative to its value in p-p or 
p-A collisions |^. Thermal parameters T = 145 MeV and fis = 370 MeV were used here and 
assumed to be centrality-independent. These values are expected in Pb-Pb collisions at 40 
AGeV. Figure 3b indicates that the statistical model in the C ensemble reproduces the basic 
features of the WA97 data |^ shown in Fig. (4): the enhancement pattern and enhancement 
saturation for large Apart- Figure 4 also demonstrates a different Apart dependence of strange 
and multistrange baryons as well as a much larger enhancement at 40 A GeV than seen in WA97 
data at top SPS. The basic predictions of the canonical statistical model is that strangeness 
enhancement from p-p to A-A collisions should increase with decreasing energy. This result 
is in contrast with UrQMD finding |^ and with previous qualitative predictions that the 
strangeness enhancement is being entirely due to quark-gluon plasma formation [|I^ . 

The quantitative comparison of the canonical model with the experimental data of WA97 
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Figure 9: Transverse momentum distribution of pion, proton and kaon obtained by STAR at 



= 130 GeV 1^ 



has been discussed in The most recent resuhs of NA57 showing an abrupt change of 
the enhancement for H as seen in Fig. 4 are, however, very unhkely to be reproducible in terms 
of the canonical approach. 



5 Particle yields and energy dependence 

In the last section arguments were presented in favour of the need for that a more general 
treatment of U(l) charge conservation, based on the canonical ensemble, if one compares the 
statistical model with experimental data for particle yields in central A-A collisions at SIS 
energies or even in peripheral collisions at the SPS and RHIC A detailed analysis of the 
experimental data in heavy ion collisions from SIS through AGS has shown that the canonical 
statistical model reproduces most of the measured hadron yields. Figure 4b shows an example 
of the recent systematic study of the comparison of the statistical model with a fully integrated 
particle multiplicities data in central Au-Au and Pb-Pb collisions at beam momenta of 1.7 
A/GeV, 11.6 A/GeV (Au-AU) and 158 A/GeV (Pb-PB) The overall agreement is seen to 
be very good. 

Figure 5 shows the compilation of chemical freeze-out parameters required to reproduce 
measured particle yields at SIS, AGS, SPS and RHIC energies. The GSI/SIS results have 
the lowest freeze-out temperature and the highest baryon chemical potential. As the beam 
energy is increased a clear shift towards higher T and lower fis occurs. There is a common 
feature to all these points, namely that the average energy per hadron is approximately 1 GeV. 
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Figure 10: Slope parameters versus particle mass at SIS PU|, top SPS pDI and RHIC p7|. 



Chemical freeze-out in A-A collisions is thus reached when the energy per particle drops below 
1 GeV a.t all collision energies |T^. The above phenomenological freeze-out condition provides 
the relation between temperature and chemical potential at all collision energies. This together 
with the measured ratio of pion/participant shown in Fig. 6a establishes the energy dependence 
of the two independent thermal parameters, the temperature and baryon chemical potential. 
Consequently the definite predictions of particle excitation functions can be given in terms of 
this model. Figures 6-8 are showing statistical model results for different particle multiplicity 
ratios along the unified freeze-out curve in comparison with experimental data. 

The statistical model predicts that the particle/antiparticle ratio should be independent 
from centrality for all collision energies. Dynamically this is a rather surprising result as 
particles and their anti-particle are generally produced and absorbed in surrounding nuclear 
medium in different ways. Figure 7b represents the energy and centrality dependence of the 

/K~ ratio from SIS to RHIC. The statistical model predictions are seen in Fig. 6b to agree 
remarkably well with the data. The measured /tt^ ratio |Q is a very abruptly increasing 
function of collision energy between SIS up to top AGS. At higher energies it reaches a broad 
maximum between 20 AGeV - 40 AGeV and gradually decreases up to RHIC energy [ET|. In 



the microscopic transport models the increase of the kaon yield with collision energy is 
qualitatively expected as being due to a change in the production mechanism from associated 
to direct kaon emission. However, the hadronic cascade transport models do not, until now, 
provide quantitative explanation of the experimental data in the whole energy range. The 
statistical model in the C formulation, on the other hand, provides an excellent description 
of K/r[ midrapidity data in the whole energy range from SIS up to RHIC, as seen in Fig. 7a. 
The abrupt increase from SIS to AGS and broad maximum of this ratio are a consequence of 
the specific dependence of thermal parameters on collision energy and canonical strangeness 
suppression at SIS. In general, however, results with the statistical model should be compared 
with 47r-integrated yields, since strangeness does not have to be conserved in a limited portion 
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Figure 11: Left: slope parameters at the SPS for different particle species. For the compilation 
of data, see for instance [pi, Right: Chemical and thermal freeze-out curve from SIS to 



RHIC. Thermal freeze-out points are from |20, 60 



of phase space. A drop in the K'^ j-n^ ratio for 47r yields has been reported from preliminary 
results of the NA49 Collaboration at 158 AGeV [0] (see Fig. 7b). This decrease is, however, not 
reproduced by the present statistical model without further modifications, e.g. by introducing 
an additional parameter 7^ ~ 0.7 |]T^ or by formi.ilating a statistical model 0/ the early stage 



m 



The appearance of the maximum in the relative strange/non-strange particle multiplicity 
ratios already seen in K^/n'^ is even more pronounced for strange baryons/meson. Figure 8 
shows the energy dependence of A/vr"*" and S~/7r"'". There is a very clear pronounced maximum 
especially in the A/vr"*". This maximum is related with a rather strong decrease of chemical 
potential coupled with an only moderate increase in associated temperature with increasing 
energy. The relative enhancement of A is stronger than that of S~. There is also a shift of 
the maximum to higher energies for particles with increasing strangeness quantum number. 
This is because the enhanced strangeness content of hadron suppresses the dependence of the 
corresponding ratio on /i^. The actual experimental data both for A/tt"'" and S~/7r+ ratios 
shown in Fig. 8 are following the predictions of the statistical model. However, as in the case 
of kaons, midrapidity results (see Fig. 8b are better reproduced by the model than Ati data (see 
Fig. 8a). 



6 Particle spectra and energy dependence 

The agreement of the equilibrium statistical model with experimental data suggests that the 
collision fireball is of thermal nature. However, this fireball is not a static object. At chemi- 
cal freeze-out the density inside a fireball is still so high that its constituents undergo strong 
rescatterings. These rescatterings result in thermodynamic pressure, which causes collective ex- 
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pansion of the medium. Consequently the particle transverse momentum distribution is showing 
a slope, which is increasing with the rest mass of the particle. This is a typical behaviour ex- 
pected from a system which undergoes transverse collective expansion. As an illustration, Fig. 9 
shows the recent RHIC results of the STAR Collaboration for the pt distribution of pion, kaon 



and proton |j57 
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Figure 12: Transverse momentum 
distribution of antiprotons obtained at 
RHIC by the PHENIX Collaboration 
for different centrality |^ . 



Figure 13: Elliptic flow parameter V2 of 
pions and protons+antiprotons in minimum 
bias collisions at RHIC |6^. Lines rep- 
resent hydrodynamical model calculations 
from 16 



The collective transverse flow is, however, not only seen in high energy data. In Fig. 10 
the dependence of the slope parameters of pion, kaon and proton is seen to be an increasing 
function of particle mass from SIS through SPS up to RHIC energy. The measured particle 
spectra at all energies seem to be affected by both thermal and transverse collective motion. If 
all particles kinetically decoupled approximately at the same time, then hadronic rrit spectra 
could be characterized by only two parameters: average thermal freeze-out temperature (Tf) 
and average flow velocity (vt) |rT|. A detailed analysis has shown that this works indeed within 
currently available data and at all energies with the exception of Q and J/ip at the SPS which 
are showing slopes lower than expected from flow systematics (see Fig. 11a. These heavy 
particles are most likely decoupled from the system already at chemical freeze-out because of 
their very low rescattering cross section with surrounding hadrons. The compilation of the 
thermal freeze-out parameters extracted from a detailed analysis of experimental data at SIS 
], the SPS and most recently at RHIC |BD] is shown in Fig. lib. At 



the AGS 



SIS, the thermal and chemical freeze-out coincide. At higher energies the thermal freeze- 
out follows the chemical one. The freeze-out temperature is remarkably consistent at AGS, 
SPS and RHIC between 100-120 MeV. Hydrodynamic evolution models pT] , |62| , |63| , however, 
reported recently that at RHIC, (Tf) ^ 128 — 140 could thus be slightly higher than at the 
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SPS. Particle pt spectra at RHIC are very well described within the concept of thermalization 
and hydrodynamical evolution. Figure 12 shows an excellent agreement of the hydrodynamical 



model 1^ with preliminary PHENIX data ||6J] for proton pt spectra measured at different 
centrality. 

The transverse flow also affects other observables, which are sensitive probes of collective mo- 
tion and equilibration in heavy ion collisions. One of these observables is the elliptic flow. This 
is particularly sensitive to the initial conditions and initial rescattering between constituents. It 
describes the azimuthal momentum-space anisotropy of particle emission in non-central heavy 
ion collisions. This emission is measured in the plane transverse to the beam direction. The 
parameter that characterizes the elliptic flow is the second coefficient V2 of an azimuthal Fourier 
decomposition of the particle pt spectra. Figure 13 shows that minimum bias data at RHIC for 



V2 are very well described by a hydrodynamical model ||6^, up to pf ~ 2 GeV. Deviations are 
seen both for higher transverse momentum and for detailed centrality analysis of experimental 
data An excellent agreement of RHIC experimental data both for particle yields and 

soft Pt spectra indicates an early thermalization of the QCD medium already on the partonic 
level and its further hydrodynamical evolution. 

7 Summary and conclusions 

We have reviewed some properties of the experimental data on particle yields and their spectra 
measured in A-A collisions from SIS up to RHIC energy. We have shown that at all energies 
particles seem to be produced according to the principle of maximal entropy, showing the 
statistical order of the multiplicities. Particle spectra, on the other hand, can be satisfactorily 
described by introducing in addition to thermal also transverse collective motion. A large 
degree of thermalization and collectivity in experimental data is particularly evident at RHIC 
and the SPS. Here chemical freeze-out conditions are remarkably consistent with those expected 
for deconfinement, and particle spectra are well described by transverse collective flow. Until 
now, however, there is no rigorous theoretical understanding of the thermal nature of particle 
production in heavy ion collisions. At RHIC and the SPS the appearance of the QGP in 
the initial state could be a driving force towards equilibration. At low collision energies the 
necessary conditions for deconfinement are not satisfied; thermalization is thus most likely to 
take place through production and rescattering of hadronic constituents. 
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